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Abstract: Improving the energy-efficiency of heating, ventilation, and air-conditioning (HVAC) 
systems has tlic potential to realize large economic and societal bcnic^fits. This paper concerns 
the system identification of a hybrid system model of a building-wide HVAC system and 
its subsequent control using a hybrid system formulation of learning-based model predictive 
control (LBMPC). Here, the learning refers to model updates to the hybrid system model that 
incorporate the heating effects due to occupancy, solar effects, outside air temperature (OAT), 
and equipment, in addition to integrator dynamics inherently present in low-level control. 
Though we make significant modeling simplifications, our corresponding controller that uses 
this model is able to experimentally achieve a large reduction in energy usage without any 
degradations in occupant comfort. It is in this way that we justify the modeling simplifications 
that we have made. Wc conclude by presenting results from experiments on our building HVAC 
testbed, which show an average of 1.5MWh of energy savings per day (p = 0.002) with a 95% 
confidence interval of l.OMWh to 2.1MWh of energy savings. 

Keywords: Model-based control; adaptive control; methodology; evaluation. 



1. INTRODUCTION 

Nearly 10% of greenhouse gas emissions and 25% of the 
electricity used in the United States is due to heating, 
ventilation, and air-conditioning (HVAC) systems in build- 
ings (U.S. DoE, 2009; McQuade, 2009). This has driven 
research into better control methods (e.g., (Nghiem and 
Pappas, 2011; Kelman and Borrelli, 2011; Oldewurtel 
et al., 2012; Liao et al., 2012; Aswani et al., 2011b; Ma 
et al., 2012)) that can help mitigate the negative exter- 
nalities due to the large energy consumption of HVAC, 
while still ensuring the comfort of building occupants. 
But the heterogeneity of HVAC equipment with respect 
to their physical modalities makes it difficult to develop a 
control design methodology that scales to many types of 
equipment. 

Even when new HVAC controllers are designed, experi- 
mentally comparing their efficiency and comfort in relation 
to existing controllers is difficult because of the large 
temporal variability in weather and occupancy conditions. 
Identifying energy models of HVAC equipment can be 
difficult because some equipment is designed to operate 
most efficiently at certain temperatures or settings, which 
requires extensive measurement to characterize. Moreover, 
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not all buildings have equipment to directly measure the 
energy consumption of only HVAC equipment. 



1.1 Experiments with HVAC Controllers 

Siroky et al. (2011) showed that a linear model predictive 
controller (MPC) could provide a 10-20% reduction in en- 
ergy usage of a ceiling radiant heating system, as compared 
to the default, manufacturer-provided controller. The tem- 
perature difference between a water set point and the 
outside air temperature (OAT) was used to approximate 
the energy usage of each controller. 

Aswani et al. (2011b) designed a new controller for an 
air-conditioner on a single-room testbed, which achieved 
up to 30% savings on warm days and up to 70% savings 
on cool days. Mathematical models of the temperature 
dynamics of different control schemes and their energy 
characteristics were constructed in order to allow com- 
parisons between different controllers of experimentally 
measured HVAC energy usage to simulations over identical 
weather and occupancy conditions. 

A 20% improvement in performance of thermal storage for 
campus-wide building cooling was achieved by Ma ct al. 
(2012) using better control methods. Direct energy mea- 
surements of the equipment along with a regression model 
of baseline performance were used to compare controllers. 
This analysis approach requires direct measurements, be- 
cause it can only statistically distinguish large differences. 




Fig. 1. Sutardja Dai Hall is 141, 000-square- foot building, 
and is part of our BRITE-S testbed. 

1.2 Overview 

This paper describes our design methodology for an 
energy-efficient controller of a building- wide HVAC system 
that is able to maintain comfortable occupant conditions. 
We begin by describing our HVAC testbed, before explain- 
ing the modeling procedure (cf. (Aswani et al., 2012b)) 
that was used to identify the thermal dynamics of the 
building and the HVAC system. Next, we describe a hybrid 
system (Tomlin et al., 2000) version of learning-based 
model predictive control (LBMPC) (Aswani et al., 2011a). 

LBMPC is a robust form of adaptive MPC. Compared 
to linear parameter- varying MPC (Kothare et al., 2000; 
Falcone et al., 2008), LBMPC differs in that it provides 
robustness to model changes using tube MPC (e.g., (Chisci 
et al., 2001)). Furthermore, the robust, adaptive MPC 
in (Fukushima et al., 2007; Adetola and Guay, 2011) 
use an adaptive model with an uncertainty measure to 
ensure robustness, while LBMPC uses an adaptive model 
to improve performance and a nominal model with an 
uncertainty measure to provide robustness. 

We conclude by experimentally comparing, on a building- 
wide HVAC system, our hybrid system LBMPC controller 
to the default controller. This is done using the comparison 
methodology described in (Aswani et al., 2012a), which 
uses nonparametric methods that can compute and com- 
pare quantitative metrics of energy usage and occupant 
comfort for different HVAC controllers. 

2. BRITE-S TESTBED 

The Berkeley Retrofitted and Inexpensive HVAC Testbed 
for Energy Efficiency in Sutardja Dai Hall (BRITE-S) 
platform (Krioukov et al., 2011; Aswani et al., 2012b) is 
a building-wide HVAC system that maintains the indoor 
environment of a 141,000-square-foot building, shown in 
Fig. 1, that is divided between a four-fioor nanofabrication 
laboratory (NanoLab) and seven floors of general space 
(including office space, classrooms, and a coffee shop). 
The building automation equipment can be measured and 
actuated through a BACnet protocol interface. 

The HVAC system uses a 650-ton chiller to cool water. Air- 
handler units (AHUs) with variable-frequency drive fans 
distribute air cooled by the water to variable air volume 
(VAV) boxes throughout the building. Since the NanoLab 



must operate within tight tolerances, our control design 
can only modify the operation of the general space AHUs 
and VAV boxes, with no modification of chiller settings 
that are shared between the NanoLab and general space. 

The default, manufacturer-provided controller in BRITE- 
S uses PID loops to actuate the VAV boxes and keeps a 
constant supply air temperature (SAT) within the AHUs. 
Conventional SAT reset control is not possible because sev- 
eral VAV boxes provide maximum air ow rates throughout 
the day for nearly the entire range of SATs. These zones 
are often dominated by heating from computer equipment. 

3. IDENTIFYING THERMAL DYNAMICS 

For the purpose of modeling and control, we will assume 
that data is sampled every 15 minutes; angle brackets (i.e., 
(•)) denote measurements sampled at this rate. Let Ts{k) 
and To{k) be the SAT and OAT, respectively, at time k. 
Similarly, Tj (k) is the temperature in the j-th zone of 
the building at time k, for j — 1, . . . , Z zones. The VAV 
box in the j-th zone controls the zone temperature Tj by 
modulating the amount of cool air sent to the zone Fj (fc) 
and the amount that the air is reheated Rj{k). 

In general, the thermal dynamics of each zone are 

T,{k + l) = f,(Ti{k),....Tz{k),F,{k),...,Fz{k), 

Ri{k),...,Rz{k),To{k),Ts{k),0,Xy (1) 

where fj{-) is some unknown nonlinear function, O are 
variables related to occupancy, and X are variables related 
to other effects like the use of equipment, solar heating, 
etc. Some simplifying assumptions are typically made by 
(a) considering the physical adjacency of different zones 
(Kelman and Borrelli, 2011; Oldewurtel et al., 2012; Liao 
et al., 2012), and (b) assuming that the effect of occupancy 
and other factors enters additively into the dynamics 
(Aswani et al., 2011b, 2012b). After these assumptions, 
the model is 

T, (fc + 1) = (t„, (fc) , . . . , r„, (fc) , Fj (k) , (k) , 

To{k),Ts{k)^+0 + X, (2) 

where {ni, . . . , Uq} is the set of zones adjacent to j. 

Additional assumptions allow further modeling simplifica- 
tions. We assume that the zone temperatures do not vary 
significantly throughout the day, since their temperature 
is in principle being controlled by the HVAC system. 
Furthermore, we assume that the additive influence of 
the occupancy and other effects can be modeled by a 
single term qj{k). Even after making these assumptions, 
the model to be identified is nonlinear since the SAT Tg 
affects the thermal dynamics in a bilinear form (Kelman 
and BorreUi, 2011; Oldewurtel et al., 2012). 

We take a hybrid system approach by forcing the SATs 
to belong to a finite set of values M — {Tg-^, . . . ,Tsj,}, 
where p is the number of modes (Aswani et al., 2012b). 
This allows us to consider multiple linearizations of (2). In 
our application to the BRITE-S testbed, we took p = 3 
and M = {52°F, 58°F, 62°F}. For the m,-th mode (for 
nii G {1, . . . ,p}) of fixed SAT T^^^ , the model is given by 



T,-(fc + l)=a™;T„,(fc) + ...a™;r„,(fc) 

+ F, (k) + cY' R, (k) + d™- To {k) + (k) , (3) 

where the coefficients are unknown scalars and qj is an 
unknown function of time k. The purpose of the system 
identification is to compute these unknown values. 



3. 1 Experiments for System Identification 

In (Aswani et al., 2012b), we used one week of data with a 
fixed SAT to identify a model. Identifying a hybrid system 
model where the SAT is able to change is more challenging, 
because identifying a model of equal fidelity would require 
three weeks (since we have three modes) of experimental 
data. As a result, we modified our modeling approach. We 
conducted experiments in which a small amount of data 
was gathered by cycling through all of the SATs in M so 
as to cover all of the modes of our hybrid system. 

Our experiment was as follows: Starting at midnight, wc 
set the SAT to T^, = 52°F for two hours. We next set 
the SAT to = 58°F for two hours. After this, the SAT 
was changed to T^^ = 64°F for two hours. During these six, 
consecutive hours, the other HVAC configuration was kept 
fixed. The reason for picking a relatively quick horizon for 
all three experiments is that the heating load is roughly 
constant over a short time span. 



3.2 Initial Parameter Identification 

We used a small amount of training data to construct an 
approximate initial model that was used to do control. To 
improve the controller performance, we have the option to 
re-identify the model using the semiparametric regression 
approach discussed in (Aswani et al., 2012b). 

The approximate initial model was constructed as follows: 
We begin by making additional modeling simplifications. 
Specifically, the exogenous heating load term qj (k) was 
changed to also include the effects of OAT and adjacent 
zone temperatures. For modeling purposes, we assume that 
the heating load term docs not change significantly over 
a short time. This was ensured by conducting the experi- 
ments for modeling within a quick time span. We know a 
priori that this last assumption is only approximate, but 
it serves to provide an initial model for which additional 
measurements can then be used to improve it. 

Suppose wc have (a) measurements for the mi-th mode 
at times k such that < k < Urm', and (b) prior 

distributions for the coefficients a™* ^/'(a™*, a™*), 

bj' ~ J^ibJ' , 6™^ ) , ~ Af , c™^ ) , where the notation 
Af{fi, S) denotes a jointly Gaussian random variable with 
mean jj, and covariance E. Our initial model is given by 

T,{k + l)=aJ-Tj{k)+b]^'F,{k)+c]^^R,{k)+q,{k), (4) 

and the coefficients can be identified by solving the follow- 
ing Bayesian, constrained least squares problem 



minE^,=i E,^ri:;(r,(fc + 1) - af^T,{k) (5) 
-bJ^^Fj{k)-cf^Rj{k)-qj{k)r 

+ (ap - ap)V57' + (6f - 67')V6f 

s.t. aj = aj^ Vr,s e (6) 

6r^<(T.,+,/T.J-6;:, Vre{l,...,p-1} (7) 

c5 = c^^ Vr,se{l,...,p} (8) 

qj{k)=qj{q), Vfc, g € Ure{i,...,p}[^r, f/r] (9) 

The constraints in the optimization problem reflect con- 
straints between different hybrid modes of the HVAC 
system. Constraint (6) ensures that the time-constants 
of the thermal dynamics are constant for each mode of 
the hybrid system, while constraint (7) encodes the fact 
that cooler temperatures will provide greater amounts of 
cooling in each zone. Constraint (8) reflects that the re- 
heating capability of each VAV box is relatively constant 
across different SATs. Lastly, constraint (9) represents the 
approximation that the occupancy is flxed over a short 
window of time. 

3.3 Modeling the VAV Box Control 

Each VAV box uses air flow Fj and reheat amount Rj 
to modulate the temperature of the j-th zone. The VAV 

boxes use a PID controller (note that the use of a PID 
controller in zone control is typical for building HVAC 
systems (Honeywell, 1997)), and so we approximate this 
as a proportional controller. This is not restrictive because 
the "learning" portion of our LBMPC controller can 
compensate for the unmodeled integrator portion of this 
control. We will discuss this in the next section. 

Let ej{k) = Tj{k) — Sj{k) be the difference between the 
zone temperature and the temperature set point for the 
j-th zone {Sj). Then, we use the following model for the 
control of reheat amount 

riOO, ffej(fc)<-l°F 
Rj{k) = } -100 • ej{k), if - 1°F < ej{k) < 0°F . (10) 
io, if ej(fc) > 0°F 

We assumed that each VAV box has the same controller, 

and so this controller model was identified by fitting 
a piecewise linear model to the observed data points 
{ej{k) , Rj{k)) over all zones. A comparison of the data 
points for a single zone and the fitted model is shown in 
Fig. 2. 

A similar model is used for the air flow amount 

(aj, ifej(fc)<0°F 
Fj{k) = < (Wj - aj) ■ ej{k) + a^, if 0°F < (k) < 1°F , 
[ijj, ifej(fc)>l°F 

(11) 

where aj and coj are the minimtun and maximum amount 
of air flow allowed in the j-th zone. These values are 
configured by the building manager. 

4. THE HYBRID SYSTEM LBMPC CONTROLLER 

Several inputs can be used for control: The reheat Rj and 
air flow Fj in each zone can be explicitly actuated, or they 
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Fig. 2. A scatter plot of data taken from a single zone of 
temperature error and reheat amount in percentage is 
shown, along with the model (solid line) that we use 
for the VAV box controller. The discrepancy is largely 
due to the integrator term that we leave unmodeled; 
our approach is to "learn" the value of the integrator 
term as we are doing control. 

can be implicitly actuated by varying the set point Sj. 
Here, we only actuate the SAT; our methodology can be 
extended to also utilize the other inputs. 

For the sake of argument, suppose that the sequence of 
SATs is fixed over a horizon of length iV and starting at 
time k = 1. So if the mode sequence of SATs is M = 
mi, . . . , TTT-AT, then the corresponding SATs are Ts{l) = 
Tg^^ , . . . ,Ts (N) = Tg^^ . Consequently, the temperature 
in each zone at time k can be computed by just "running" 
or simulating the model forward in time; no optimization 
is needed to compute the values of Fj,Rj,Tj, etc. 

As a result, minimizing our cost function subject to the 
thermal dynamics and system constraints requires opti- 
mization over only the set of possible sequences. In other 
words, our optimization problem is an integer program. 
Fortunately, there is a constraint that greatly reduces 
the computational complexity: The SAT can only change 
once every hour; the reason for this constraint is that 
large, frequent changes to the SAT can damage the HVAC 
equipment. So if the horizon is = 16, which corresponds 
to four hours when sampling at 15 minute intervals, then 
there are 3621 different combinations. 

In order to further simplify the computational complexity, 
we make use of a heuristic that further reduces the 
computation. We specify that the mode is fixed over the 
span of every hour. For our setup — where we have three 
modes and four possible mode changes — this means 
that we have to compute the cost for 3^ = 81 different 
combinations. This is a low number of combinations that 
can be computed under one second on a desktop computer 
because, as mentioned earlier, we do not need to optimize 
over other variables. For future extensions where other 
variables are used to do control, we would only need 
to solve 81 convex optimization problems (specifically 
quadratic programs), which can be reasonably solved. 

4-.1 Energy Modeling 

Before we present the optimization formulation of the 
controller, we provide some intuition into the form of the 
cost function we use. A building HVAC system typically 
has several individual pieces of equipment that contribute 
to the overall energy usage. Within BRITE-S, most of the 



energy consumption is due to three elements: the fans in 
the AHU, the chiller that cools water and indirectly cools 
air, and the reheating that occurs each zone's VAV box. 

Even though parameterizations of the energy usage of 
the equipment are known (Kelman and Borrelli, 2011; 
Oldewurtel et al., 2012), modeling these features is difficult 
because individual energy measurements are usually not 
available. There is another subtle point regarding the 
energy models and their relationship to our cost function. 
Let Ei{-) , E2{-) T E^{-) be the energy due to fans, chiller, 
and reheating. One possible cost is '}Zj{'^3~ ^j)"^ 
HE2{-) + 7£'3(-), where A,/x,7 G R are constants. There 
are two reasons for weighting energy usage based on type: 
First, the energy usage of the equipment is often related to 
the mechanical and physical stress on the equipment, and 
so differential weighting allows a finer level of regulation 
with respect to such considerations. Second, the energy 
usage of equipment sometimes comes from different energy 
sources. For example, heating is provided by steam while 
the fans are powered by electricity in BRITE-S. 

Because the cost function weights each energy usage func- 
tion, we really only need to know the energy model up to 
a constant, unknown scaling factor, because this gets sub- 
sumed into the scaling in the cost function. This simplifies 
the energy modeling that we do for the purpose of control. 
We use the following energy models: The fan energy usage 
is El ^ iJ2j Pj)^ ^-nd the chiller energy is E2 ~ {To — 
Ts) ■ J2j Pj (Kelman and Borrelh, 2011; Oldewurtel et al., 
2012), while the reheat energy is E^ ^ J^j Rj- 



4.2 Optimization Formulation 

We can now present the optimization formulation of the 
hybrid system LBMPC controller. The basic intuition 
behind LBMPC (Aswani et al., 2011a) is that two models 
of the system are kept. The first is a nominal model that 
is used with respect to the constraints, and the second 
model is used in the cost function and is updated using 
data gathered during control. This maintains robustness 
while improving performance through model updates. 

Without loss of generality, we assume that the control 
action is being computed for time k = 1; also recall 
that we do control at a rate of every 15 minutes. If 
Tj{l) , Rj{l) , Fj (1) are the predictions of the linear model 

from time fc = 0, then let qj{i} = T,(l) - Tj(l), fj{i) = 
Fj{l) -Fj (l), f-j{i) = i?j (l) -Rj{l). The intuition is that 
these terms with hats represent corrections to the predic- 
tions of the MFC and provide the adaptation inherent in 
the controller (Aswani et al., 2011b). Specifically, qj rep- 
resents an estimate of the heating load due to occupants, 
weather, solar heating, and equipment; fj represents the 
integrator term in the PID control of air flow in the j-th 
VAV box; and fj represents the integrator term in the PID 
control of reheat amount in the j-th VAV box. 

The control action is given by the minimizer to 



min Eli(E,m + l)-Sj{i + l)r (12) 

mi,...,mjv \ 

+ MM)-T,(i)).E,-^;w) 

s.t.(4),(10),(ll) 

Fj = Fj + fj ; Rj — Rj + fj 
f,{z + 1) = «7'T,(z) + bY'F,{^ + c7-i?,(z) 
+ qj{i) 

T,{Aq + r) =T,{4q + s), Vr, s € {1, . . . , 4} (13) 

A g G {0, . . . , riV/41 - 1} 
Tj{i) e [66°F,78°F], yjG{l,...,Z} 

Note that constraint (13) allows the SAT Tg to switch 
value only once an hour, which reduces the computational 
complexity of the controller as discussed previously. 

A desktop computer took an average of under one second 
to solve this optimization problem for BRITE-S. Further- 
more, we used values of A = 6.7e4/(^^ ctj)^, M = 1.3e— 3, 
and 7 = 6.7. These values were picked using an iterative 
process in which (a) the control was computed but not 
used for actuation; (b) the control and predictions were 
analyzed for if SAT stayed at higher values, the reheat 
amount was low, and the air flow remained at moderate 
levels; (c) the coefficients A, /x,7 were changed and the 
process starting at (a) was repeated until (b) was satisfied. 

5. MEASURING EFFICIENCY 

Separate energy measurements of the general space HVAC 
system are not available in BRITE-S. Instead, we only 
have access to measurements of both the general space 
HVAC and the NanoLab HVAC. We denote these energy 
measurements as E[i] for i = 1,...,D, where the index 
i is over hourly intervals. Furthermore, we have measure- 
ments of the OAT that correspond to the energy usage 
measurements: To[i] ior i = 1, . . . , D. 

The general model describing the relationship between 
energy usage, OAT, occupancy O, and other factors (e.g., 
solar effects, equipment, etc.) X is E = f{To, 0,X), where 
/(•,-,•) is an unknown, nonlinear relationship. But because 
occupancy and other factors are generally not directly 
measured, it is typically only possible to consider the 
energy usage averaged over occiipancy and other factors 
EoA^o) = V.[f{To,0,X)\To]. This can be estimated 
using nonparametric regression (Gyorfi, 2002), and it is a 
curve that describes the relationship between the OAT and 
the average energy consumption. Intuitively, if the data 
points (To[i], E[i\) for i = 1, . . . , D represent a scatter plot 
of energy usage versus OAT; then Eo^^iTo) represents the 
smoothed version of the scatter plot. 

The average amount of energy used in one hour is therefore 
Eg = Eg{EoATo)) = Eg [E[/(r„, O, X)\To]], where g{To) 
is a probability distribution of OATs. This notation allows 
us to define the average amount of energy used in one 
day as E^ay = J2^ti Eg., where gi{To) is the probability 
distribution of OAT during the i-th hour of the day. For 
simplicity, we assume a uniform distribution for the OAT. 



This allows us to compare the energy efficiency of two 
controllers. We can compute the quantities defined above 
for each controller, and a nonparametric methodology 
(Vilar-Fernndez et al., 2007) can be used to determine 
whether the differences are statistically significant. 

6. MEASURING OCCUPANT COMFORT 

In order to define a quantification of comfort that is both 
tractable and will scale to many buildings, we focus on 
a measure that is only dependent on temperature. We 
assume that the average temperature for the j-th zone 
is measured at hourly intervals Tj[i] iov i = 1, . . . , D. 

Let (.x)_|- be the thresholding function, which is defined 
so that (a;)-|- = if a; < and {x)+ = x otherwise. We 
define our quantification of comfort using soft thresholding 
as C = 1/Z ■ J2f=iIoi\Tj - Sj\ - Bj)+dt, where the 
integral with respect to dt is over one hour of time, 5^- 
is the set point of the zone, and Bj is the amount of 
temperature deviation for which the building has been 
configured. The intuition is that this quantity increases 
whenever Tj exceeds Sj by more than Bj , and the amount 
of increase in this quantity is proportional to the amount 
and duration of temperature deviation. The BRITE-S 
building is configured for Bj = 1°F for all zones and times. 

7. EXPERIMENTAL RESULTS 

The hybrid system LBMPC controller was used to control 
the SAT in BRITE-S for 8 days that spanned both 
weekdays and weekends, and this was compared to 22 days 
in which the default, manufacturer-provided controller 
was used. Recall that our comparison methods implicitly 
account for variations in energy usage due to occupancy, 
and this is made explicit in the 95% confidence interval for 
estimated values. The building and HVAC configurations 
were kept identical for when both controllers were used, 
and this configuration has been in use for over one year. 

For notational reasons, we use a superscript 1 to refer 
to the default controller, and superscript 2 denotes the 
hybrid system LBMPC controller. Also, we define the 
quantity AE'^'^ = E'^^^ — E\^y to be the estimated 
difference between the average energy usage over a day of 
the LBMPC controller and that of the default controller. 
The value A(7^'^ = C!day ~ ^day is an analogous quantity 
of the estimated difference in average comfort over a day. 

The estimated energy characteristics El ^{Tq), E'^ ^(Tq) 
are shown in Fig. 3, and their differences are statistically 
significant {p = 0.002). The estimated difference in average 
energy usage over one day AE'^'^ = — 1.53MWh is statis- 
tically significant {p = 0.002). And the 95% confidence 
interval is A£;2,i g [-2.07, -1.02]MWh. 

The estimated comfort characteristics (^^ ^.(To),^^ ^.(To) 
arc shown in Fig. 4, and their differences arc not statisti- 
cally significant (p = 0.8). Furthermore, the estimated dif- 
ference in average comfort over one day AC'^'^ = — 0.75°F 
is not statistically significant (p = 0.5), meaning that there 
is not enough evidence to exclude that AC^'^ = 0°F. 

The LBMPC controller provides modest energy savings at 
most OATs, which sum up to significant savings over a day. 
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Fig. 3. The cross marks and solid line (points and dashed 
line) denote the energy characteristics of the LBMPC 
(default) controller. 
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Fig. 4. The cross marks and solid line (points and 
dashed line) denote the comfort characteristics of the 
LBMPC (default) controller. 

And because the difference in comfort characteristics is not 
statistically significant, this suggests that the LBMPC and 
default controllers provides comparable levels of comfort. 

8. CONCLUSION 

We have presented a hybrid model of building HVAC and 
described its control using hybrid system LBMPC. Ex- 
periments show substantial savings, and future directions 
for further energy savings were discussed. More broadly 
speaking, our experiments on BRITE-S, and previously on 
BRITE, show that the LBMPC methodology can provide 
significant energy savings for a wide variety of HVAC 
systems operating using different physical modalities. 
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